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Thermodynamic and dynamical properties of filling-control metal-insulator transition (MIT) 

^js. \ in the Hubbard model are studied by the operator projection method, especially in two dimen- 

fvj , sions. This is a non-perturbative analytic approach to many-body systems. The present theory 

incorporates the Mott-Hubbard, Brinkman-Rice and Slater pictures of the MIT into a unified 

framework, together with reproducing low-energy narrow band arising from spin-charge fiuctu- 

TJ^ I ations. At half filling, single-particle spectra A{lj, k) show formation of two Hubbard bands and 

I , their antiferromagnetic shadows separated by a Mott gap in the plane of energy ui and momen- 

5-H ■ tum k with lowering temperatures. These four bands produce splitting to two low-energy narrow 

j/j ' bands and two SDW-like bands in the dispersion. Near half filling, the low-energy narrow band 

persists at low temperatures. This narrow band has a particularly weak dispersion and large 

j^ , weights around (vr, 0) and (0, n) momenta. The velocity of these low-energy excitations is shown 

to vanish towards the MIT, indicating the mass divergence as in the Brinkman-Rice picture, 

but most prominently around (vr, 0) and (0, n) with strong momentum dependence. This reflects 

FT-j ' the suppression of the coherence near the MIT. Main structures in A[uLi,k) show remarkable 

f— I , similarities to quantum Monte-Carlo results in two dimensions as well as in the one-dimensional 

Q . Hubbard model. The charge compressibility appears to diverge with decreasing doping concen- 

fj ' tration in both one and two dimensions in agreement with the exact and quantum Monte-Carlo 

results. We also discuss implications of the fiat dispersion formed near the Fermi level to the 

observations in high-Tc cuprate superconductors. 
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f^ ■ §1. Introduction 

The issue of the metal-insulator transitions (MIT) postulated by MotttP has opened extensi]i£j_studies 
on the subject from both theoretical and experimental aspects.H-' Pioneering works by HubbardLfQ'EF and 
j^ ■ KanamoriEl' played an active part in theoretical studies on strong electron correlations, in narrow-band 

electron systems.LF In spite of a long history of studies on strongly correlated electrons JU'EP an appropriate 
and satisfactory theoretical description of the MIT remains controversial, even for the simple Hubbard 
' O ' model. Particularly, in one and two dimensions, fluctuations become significant near the MIT. This makes it 

o 
o 



^ ■ difficult to reliably descrihe-physical properties related to the MIT. This paper aims at developing a formalism 

^\ \ of operator projectionErcl'113'' to construct a new non-perturbative analytic theory of many-body correlated 



electrons. Ej-* This method.allows us to systematically improve th£ Hartree-Fock theory, perturbation theories, 
Hubbard approximationaj'Lf and two-pole approximations. lj'lj' Using this new method, we give static and 
dynamical properties of the Hubbard model in one and two dimensions at and near half filling. 
JH \ Early theories on transitions from metals to insulators in clean systems are divided into three pictures: (i) 

By means of the Green's function technique based on Heisenberg equations of motionjlil) Hubbard described 
an MIT as the splitting |X)f two Hubbard bands due to strong local Coulomb repulsion compared with 
the original bandwidth.ola^ This description is based on the view provided by Mott.M Tho-Hubbard I 
approximation was improved as two-pole approximations in a spirit of operator projections.li3't3 These 
approximations are valid in the strongly correlated limits. However, they violate the Luttinger jSum rule 
in paramagnetic metals. They do not reproduce the correct low-density limit given by Kanamoriu' either, 
(ii) Using the Gutzwiller variational methodjlj) Brinkman and Rice predicted that at half filling, an MIT 
occurs at a finite critical interaction strength Uc due ta_thc disappearance of Fermi-liquid quasiparticles 
accompanied by the divergence of the quasiparticle mass.ll3' However, spin correlations as well as incoherent 
excitationSpappearing as the Hubbard bands are neglected in this theory, (iii) Early arguments on the Slater 
mechanismllJ) of transitions from paramagnetic metals to antiferromagnetic (AF) insulators are based on 



the Hartree-Fock theory. In this mechanism, the folding of the Brillouin zone accompanied by the AF phase 
transition produces an energy gap in single-particle excitations. Gaussian fluctuations around the Hartree- 
Fock solution are treated by a random phase approximation (RPA). Self-consistent renormalization (SCR) 
theory developed by Moriya ancLpoworkers modified the Hartree-Fock approximation and the RPA in the 
problem of itinerant magnetism.IlM' However, these theories fail to reproduce Mott-insulating features which 
persist even above the Neel temperature Tn . i-i r. 

Experimentally, Mott-Hubbard MITs have been observed in a lot of compounds,El'0^ e.g., V2O3, 
Lai_a;Sra;Ti03 and Ji^i. ^-Sr^VOg Especially, V203-based compounds have provided a typical example for 
studying the MIT.Ej'ErB' In a gross view, as the temperature and the bandwidth are varied, the compounds 
exhibit phase transitions among three different phases; antiferromagnetic insulating (AFI), paramagnetic 
metallic (PM) and paramagnetic insulating (PI) phases. The ground state is AFI under weak pressure and 
PM under strong pressure. The AFI phase has been considered to originate from the Slater mechanism. 
On the other hand, a strong renormalization of the quasiparticle mass occurs in the PM phase, as in the 
Brinkman-Rice picture. With increasing temperature, the AFI system undergoes an AFI-PM phase transi- 
tion under strong pressure, while a PI phase appears above Tn under weak pressure. Electronic properties in 
the PI phase are consistent with the Mptt-Hubbard picture. This compound also undergoes a filling-control 
MIT with decreasing V concentration.tB' In spite of a long history of the studies, a sufficient unified de- 
scription of overall properties in the global phase diagram is still missing. Lai_j;Srj;Ti03 and Lai_j;Srj;V03 
are other examples where a filling-control MIT occurs accompanied by the AF phase transition and the 
quasiparticle mass enhancement. These are other examples which can not be described without overcoming 
the three separate scenarios of the MIT. j— . 

The discovery of high-Tc cuprate superconductorsEiP has promoted further intensive studies on the is 
of the MIT with the aim of explaining their peculiar features observed near the Mott-insulating phascc. 
A proper understanding and description of the high-Tc cuprate superconductivity remains still challenging, 
since it is difficult to correctly describe strong electron correlations that prominently emerge near the MIT. 
High-Tc cuprate superconductors have a quasi two-dimensional electronic structure. Anderson proposed 
that their effective model is giiien by the single-band Hubbard model for the Cu 2>d^2_y2 orbital strongly 
hybridized with O 2p orbitalsBj) The undoped compounds are Mott insulators and exhibit the AF long- 
ranged order at low temperatures. Holes doped into these Mott insulators induce superconductivity at low 
temperatures. With further increasing holes or raising the temperature, they undergo a phase transition 
from the superconducting (SC) phase to a paramagnetic metal. Suppression of the quantum-mechanical 
coherence seems to reduce T^ with decreasing hole concentration below the optimum value. In contrast with 
the normal metals with a reduced Fermi temperature, the simple Brinkman-Rice picture does not describe 
the mechanism of this suppression of the rnhnrnnrn This is because a pseudogap appears in low-energy 
single-particle, spin and charge excitations.oEj'cSEZfEa'' AF spin excitations have a low-energy mode. Its 
characteristic energy decreases with decreasing T^ within the pseudogap. The competition of this mode with 
the singlet-pair correlations giving rise to the superconductivity appears to be the underlying origin of the 
pseudogap behavior. Pseudogap in charge excitations seems to continuously evolve into the SC gap which 
does not scale with T^. Towards the AFI phase, spectral weights outside the pseudogap must shift to high 
energies above the Mott gap. However, there seems no reliable theoretical framework that can describe such 
strong correlation effects near the MIT. For this purpose, it is crucial to develop a theoretical framework 
which gives a unified scheme containing the three aspects of the MIT, at least. 

Therefore, construction of a unified theory which captures all the above three pictures of the MIT has 
attracted a general interest in the field of correlated electrons not only in the context of the specific subject 
of the MIT. In fact, any one of the above three pictures alone is not sufficient to describe the MIT even in the 
Hubbard model. Without a unified scheme, one can not reproduce the global phase diagrams of compounds 
that undergo an MIT. It is also difficult to clarify the mechanism of the high-Tc cuprate superconductivity. 

The Bonkman-Rice solution is also derived as a-mmma gnetic saddle point in a functional integral repre- 
sentationcS based on the slave-boson formulationJ^JH-BESP namely slave-boson mean-field (SBMF) theory. 
Including the Gaussian fluctuations around thiSj-aonmagnetic saddle point gives the Hubbard bands and 
incoherent features of single-particle excitations.t3'E3' Short-ranged spin correlations, however, appear as 
higher-order loop corrections and were neglected in these theories. Instead, by considering both nonmag- 
netic and magnetic saddle points, the functional integral theory based qjo, the slave-boson formulation gives 
an overall treatment of ferromagnetism, antiferromagnetism and MIT.E2P The results agree with those of 
Hartree-Fock theory in the weak-coupling regime. They are also consistent with Nagaoka's theoremEj* and 



Kanamori's resultaS* on ferromagnetism in the strong-coupling regime. For further improvements to in- 
clude short-ranged spin correlations, however, it is required to overcome a difficulty associated with the 
uncontrolled treatment of the local constraints, j— . 

The dynamical mean- field theory (DMFT)E3'lZP succeeded in incorporating the Hubbard and Brinkman- 
Rice pictures into one theoretical framework. This theory gives the exact Green's function in infinite dimen- 
sions. It predicts that at half filling, a resonance appears at the chemical potential inside the Hubbard gap 
below a critical interaction strength Uc- The resonance has a constant intensity and a decreasing bandwidth 
as U increases. This agrees with the Brinkman-Rice picture.E^ The DMFT has also been applied to the 
finite-dimensional Hubbard model and other correlated models. H However, it suffers from the following two 
related oversimplifications. Though this theory gives correct dynamical properties of the local self-energy 
part, it neglects its momentum dependence. In addition, the DMFT ignores the short-ranged two-particle 
correlations including spin correlations. 

Some attempts have been made to improve the dynamical mean-field tbenry (T^MFT) Recently, Jarrell 
and coworkers have developed the dynamical cluster approximation (DCA).e3oI1J'OOcJ-' In the DCA, the 
cluster probleiiLwith A^c(> 1) sites is solved. In the case of A^c ~ 1, the DCA gives the original formulation 
of the DMFT.H Then, N^ independent momentum regions appear in calculating the self-energy part. For 
the two-dimensional Hubbard model at half filling, the DCA succeeded in reproducing a Mott insulator with 
sharp peaks at the low-energy edges of the Hubbard bands energetically separated by an antiferromagnetic 
(AF) pseudogap at finite temperatures. At weak-to-intermediate couplings, this prediction contrasts with 
the paijainagnetic solution of the DMFT. The DMFT predicts that the resonance emerges at the lowest 
energy .£3 In the strong-coupling regime, the DCA and DMFT give similar Mott-insulating features. In the 
DCA, however, one has to solve the cluster problem. Although it captures some aspect of strong momentum 
dependence of the self-energy part, this gives a severe limitation due to a finite cluster size Nc- This makes 
it difficult to consider the thermodynamic limit. It is also difficult to study more complicated systems with 
many degrees of freedom. Therefore, it is desired to build another analytic theory for studying correlated 
many-electron systems so that neither the system size problem nor the cluster size problem occurs. 

Froiiua—viewpoint of the Slater picture, by using the diagrammatic methocLand/nr. Luttiiiger.-Wfird func- 



tional,Ea'C3'CII' microscopic treatmcptaJ:o modify the RPA have been made.L§L3'l§oE3H'E§E§ES' They 
involve conserving approximations ,£3^113^ the fluctuation-exchange api)p:pdmation,E3' paramagnon theories, 
parquet approachescH' and the two-patticle self-consistent method.EEfE^P They can describ£_.Fermi-liquid 
states in the weak-correlation regipaH while they fail to predict Mott-insulating featuresJ^S^ except the 
two-particle self-consistent met hod. El' S The two-particle self-consisten|Ljn£.thod reproduces the AF pseudo- 
gap in the half-filled square-lattice Hubbard model at low temperatures.E^EiP However, the Green's function 
and the self-energy part are not self-consistently determined in this method. 

Apart from analytic theories, the quantum Monte-Carlo (QMC) simulation is a powerful tool to study 



properties of the MIT in 



strongly correlated models if the negative sign problem is not serious.B' The st 
the two-dimensional Hubbard model were investigated by QMC calculations .Ej'E3'EIP The compressibility, 
thus the charge susceptibility was shown to diverge towards the MIT.H QMC calculations, supplemented 
with an approximate analytic continuation by the m ay irn iim en tropy method, have given some characteristic 
dynamical properties of the MIT in two dimensions:E3'E3'E2HEH£3' Nearly separated Hubbard bands develop 
near half filling. The Luttinger sum rule which is a criterion of the validity of the Fermi liquid in two- and 
three-dimensional metals at T = appears to be violated at least at T > J = At^ /U possibly due to thermal 
effects.H This reflects that the Fermi temperature is suppressed near the MIT. Near half filUftg, a sharp peak 
of the local density of states develops near the Fermi level with decreasing temperaturcES Furthermore, 
at low temperatures T < J, extra two bands appear in addition to the Hubbard bands-, A strongly weak 
dispersion around the (tt, 0) and (0, tt) momenta has been found at zero tercifje£ature.EJ These features 
inferred from QMC calculations agree with the results of exact diagonalization.Ej'Ej) The strong momentum 
dependence of quasiparticlc excitations is a crucial feature to be reproduced beyond the Mott-Hubbard, 
Brinkman-Rice and Slater pictures. 

We next consider oneidimeiisional cases. Usual one-dimensional metals are classified into the Tomonaga- 
Luttinger (T-L) liquidJ£j'E3LZP Spinous and holons play a role of elementary excitations in the T-L liquid. 
The Bethe ansatz wave functionO for the Hubbard model in one dimension suggests that in the case of 
[/ = oo, spinops and holons are completely decoupled, and thereby the spin-charge separation is realized 
at all energies.EB' However, available spectral properties of the one-dimensional Hubbard model that can 
be extracted from the wave function are severely restricted to the case of half fiUingJlj) one hole doped 



into thc-Mptt insulatorjO) the strong-coupling limit ,00 the close vicinity of the Mott insulatorO or 
uj = Ox3l3' For finite values of U, a mixing of spin and charge collective degrees of freedom occurs. In 
this case, information involved in the Bethe ansatz wave functionH is not sufRcient for obtaining one- 
electron excitation spectra. Instead, QMC calculations supplemented with the maximum entropy methodic' 
have given onctclectron spectral properties at half filling which agree with those obtained by the SBMF 
approximationHj' For the one-dimensional t-J model, exacLxliagonalization studies indicate the presence 
of spiniin_and holon branches in one-electron excitationalU'O that are expected from the arguments for 
J = OSlIj) One-electron excitations are formed by the convolution of spinous and holons. The bandwidth 
is given by ttJ for spinous and At for holons. Then, they compose two branches in one-electron excitations 
mainly characterized by spinous and holons. At half filling, the holons become massive and this produces a 
Mott gap in the one-electron excitations. When the spin-charge separation occurs prominently, the spinon 
branch in the lower Hubbard band at half filling is found only for |fc| < 7r/2. Then, the complete shadow 
structure is missing in one-electron excitations. This is allowed in one dimension, because even at T = 0, there 
is no AF long-ranged order and the translational symmetry is never broken. For the Hubbard model with a 
finite value of U, however, a coupling between spinous and holons makes it difficult to discuss properties of 
one-electron excitation spectra. 

In spite of these studies mentioned above, a sufficient theoretical framework that gives a unified treatment 
of Hubbard, Brinkman-Rice and Slater pictures together with strong momentum dependence of quasiparti- 
cle excitations and reproduces key features obtained in numerical simulations is still missing. Particularly, 
it remains unsolved how electronic spectra change from metals to Mott insulators. This is an important 
subject in the investigation of the MIT. To study this subject, in the rest of this section, we introduce a 
promising theoretical formulation of the operator projection method (0PM) employed in the previoustiP and 
the present papers. The 0PM is a non-pcrturbative analytic approach which systematically improves the 
Hartr-efi-pEhcV. theniy. ±he Dyson-equation formalism. of the Green's function, theories-based on spin fluctua- 
tionsJl3'E3'C§EJ'E3'EilE3' Hubbard approximationaS'la) and two-pole approximations J13'0E2P The 0PM also 
reproduces a sharpj)eak near the Fermi level in doped cases, in agreement with the results of the dynamical 
mean-field theory .H The velocity that characterizes the low-energy single-particle excitations is shown to 
vanish towards the MIT particularly around (tt, 0) and (0,7r). This suggests that the quasiparticle mass 
diverges towards the MIT, though the enhancement of the damping rate makes it difficult to define quasi- 
particles. This mass divergence is to some extent consistent with the Brinkman-Rice picture, although its 
strong momentum dependence is beyond the Brinkman-Rice picture. Moreover, the 0PM reproduces a sharp 
peak of the density of states at each edge of the Mott gap that leads to the divergence of the compressibility 
towards half filling. This contrasts with the Brinkman-Rice and dynamical mean-field theory. Therefore, 
the method opens a way to construct a unified scheme beyond the three different pictures. Furthermore, our 
theory given in the previoustj-' and the present papers does not suffer from the problem on the system size, 
unlike the dynamical cluster approximation.cS' 

The|43rojection technique in the operator space to calculata correlation functions in the linear response 
theorylE|) was developed by NakajiroajB ZwanzigEP and MoriO This is a non-perturbative approach based 
on Heisenberg equations of motiono^ It leads to a continued-fraction expansion of a correlation function. 
The expansion is uniquely determined and can also be described by the series of Dyson equations. This 
continued-fraction expression also uniquely determines a corresponding moment expansion. Though an 
exact proof is absent, the continued- fraction expansion is empirically known to converge more rapidly than 
series expansions and to have a larger convergence region. Any theory involving an equation of motion is 
called an equation-of-motion approach. _Jn_i|ius_scH.sc, the operator projection theory belongs to a family 
of equation-of-motion approachcs.EHa'tl't3'E3E]'E3'L3Eil' Operator projection theories are, however, more 
sophisticated in that it produces a systematic way of performing a continued-fraction expansion in energy. 
The 0PM opens a possibility to systematically extract the lower-energy spectral properties projecting out 
higher-energy degrees of freedom in a self-consistent fashion. 

The formulation of the 0PM is explicitly given in § g. Here, we summarize the connection of the 0PM 
with other known theoretical frameworks. 

If one performs the projection procedure for the single-particle operators up to the first order, it gives 
the Hartree-Fock theory, in the symmetry-broken phases as well as symmetry-unbroken phases, the free- 
electron properties. At this stage, the higher-order dynamics remains as the self-energy part in the so-called 
Dyson equation for the Greenls. function. This gives a starting point of the perturbation theory, and thus 
conserving approximations ,OeZP the fluctuation-exchange approximation Jls' paramagnon theories, parquet 



approachesc2l' and the two-particle self-consistent method.BEJ As we have already mentioned, these theories 
fail to predict Mott-insulating features in a self-consistent manner. Here, we note that this is because only the 
zeroth- and first-order moments of the Green's function are exact at this stage. To study strong correlation 
effects that prominently emerge near the MIT, a natural step is to find a way of treating the self-energy 
part in a systematic fashion beyond the Gaussian approximations. Since-the perturbation theory works for 
the self-energy part more effectively than for the Green's function itselfjlj) it seems natural to expect that 
considering the Dyson equation for the self-energy part improves its estimate. The projection technique is 
useful for this purpose. 

The second-order projection filters out the local nature of the Hubbard model—Here, the atomic properties 
and the superexchange term are involved. In the Hubbard approximations, a u'Q^ this procedure was partly 
taken: Hubbard obtained the exact second-order moments of the Green's function but not the third-order 
moments. His attempts of making a self-consistent theory do not constitute the systematic projection 
proceduifi. Roth improved the Hubbard I approximation along the concept of the operator projection 
theory.E3' This introduces the superexchange term into the self-energy part which gives rise to a correct 
treatment of the third-order moments. However, the higher-order dynamics appearing in the Dyson equation 
of the self-energy part, which is hereafter called the second-order self-energy part, was ignored. Then, his 
results for the Green's function have only two poles (two-pole approximation). Some attempts to take 
account of higher-order dynamics were made by Matsumoto and Mancini.H' In their method called the 
composite operator method, the higher-order dynamics were included only within the two-site level in a 
local picture. However, it is not known how reliable the composite operator method is in the presence of 
short-ranged fluctuations. It is crucial to incorporate short-ranged AF spin correlations near continuous 
magnetic phase transitions. These short-ranged correlations can be included in the operator projection 
method (0PM) by several choices of decoupling approximations for evaluating the second-order self-energy 
part through convolutions of the single-particle Green's function and two-particle susceptibilities. M-* 

The results obtained with the present method are given in § |3[ The present 0PM reproduces the following 
single-particle properties in two dimensions: (I) The filling-control MIT takes place at half filling. A Mott 
gap opens and completely separates the two Hubbard bands. (II) At low temperatures, strong AF spin 
fluctuations produce an AF shadow in each Hubbard band, yielding a four-band structure. This structure 
in energy and momentum is regarded as a superposition of two low-energy narrow bands and two SDW-like 
bands. (Ill) Upon hole (electron) doping, the narrow bands appear to naturally merge into the top of the 
lower (upper) Hubbard band. Then, the low-energy band persists near half filling. This narrow band has 
a particularly weak dispersion and dominant weights in the momentum regions around (tt, 0) and (0,7r). 
(IV) Due to these momentum regions, a sharp peak of the local density of states grows at the Fermi level. 
With increasing temperature, the sharp peak diminishes due to thermal effects. (V) The Fermi velocity 
characterizing the low-energy narrow band is found to vanish as (n) —^ 1, reflecting a suppressed Fermi 
degeneracy near the MIT. This suggests that the quasiparticle mass diverges towards the MIT. This is 
consistent with the Brinkman-Rice picture. However, our results show that an enhancement of the damping 
rate of quasiparticles makes it difficult to define quasiparticles. It also shows characteristic and strong 
momentum dependence of quasiparticles excitations beyond the Brinkman-Rice picture. Such enhanced 
momentum dependence may generate momentum-dependent electron differentiations in the low-energy level 
of the energy hierarchy formation. (VI) For the one-dimensional half-filled Hubbard model, the 0PM also 
reproduces a Mott insulator at an intermediate coupling. The Hubbard bands separated by a Mott gap 
appear. At lower temperatures, their AF shadows develop similarly to the two-dimensional case, in spite 
of the absence of the AF long-ranged order. The single-particle excitations in one dimension consist of the 
convolutions of the spinon and holon excitations which couple each other in the case of finite values of U . 
The present results suggest that at an intermediate coupling, the mixing of spinons and holons is too large 
to interpret the singlp-particle excitations in terms of the simple picture based on the completely. drjrouj) led 
spinons and holons.L2) These features obtained by the 0PM are consistent with QMC resultsJ^'t3Ljl£iJL3'L3' 
as we will discuss in this paper. 

The Mott insulators are incompressible, while the metals are compressible. This is reflected in the charge 
susceptibility; in the Mott insulator, the charge susceptibility vanishes. On the other hand, approaching tlue 
MIT in the metallic phase, the compressibility, thus the charge susceptibility diverges towards the MITE' 
The present method also correctly reproduces this remarkable property. 

It is remarkable that the present theory captures the three essential features of the MIT as given in 
the (i) Mott-Hubbard, (ii) Brinkman-Rice, and (iii) Slater pictures. In fact, this theory reproduces the 



Mott-Hubbard picture at moderate temperatures near the MIT. It also reproduces the enhancements of the 
quasiparticle mass and damping rate near the MIT. This reflects the suppression of the Fermi degeneracy 
near the MIT. This reminds us of the Brinkman-Rice picture. The growth of AF spin fluctuations at low 
temperatures yields AF shadows of the Hubbard bands in the single-particle dispersions. In § H, we will give 
our view of the MIT in the two-dimensional Hubbard model. All the three aspects are crucial to correctly 
describe the MIT. We also note that previous equation-of-motion approaches have not taken account of 
short-ranged AF spin correlations. 

In § pi we compare the present results with experimental observations in high-Tc cuprate superconductors. 
Similarity of the narrow band obtained here to that observed in high-Tc cuprate superconductors implies 
that to discuss the superconductivity, it is required to include strong electron correlations related to the 
MIT. 

§ is devoted to the summary of this paper. We also discuss problems left for future study. 

§2. Formulation of operator projection method 

Below, we follow the general formalism of the 0PM which has been briefly discussed in the previous 
letterO We take the Hubbard Hamiltonian with an electron transfer tx,x' from an atomic site x' to x and 
the local Coulomb repulsion U; 

H = - ^ tx,x'ci^Cx's + U^nx^nxi, (1) 

X,X',S X 

with Uxs = CxsCxs- For simplicity, we restrict ourselves to the case of tx^x' = tx',x- Applications to more 
general situations are straightforward. We define H = H — jjl 'Ylixs ^^s with the chemical potential fx. We 
perform the projection procedure for the electron creation and annihilation operators at an atomic site x 
with a spin index s, Cx^s and Cx,s- It is useful to define the fermionic 4A^-component vector operator ^, 
which is composed of 

,T, - f P0 + P3 -PO- PZ \t( t t \ 

"""i 2 "^""^ 2 '^V (ca;T'Cxi,c;^pcJ^J 

= *(Ca;T,Ca:i,-4^,4|). (2) 

in the real-space representation. N represents the number of the atoms. ^^ is defined as the Hcrmitian 
conjugate operator of ^ . In the real-space representation, it is written as 

*i == (4T'4i'-Ca;i,CccT)- (3) 

(To (po) and cr (p) are the identity and Pauli matrices, respectively, which operate to the spin (charge) space 
of vector operators as ^ from the left-hand side and those as ^^ from the right-hand side. 

We define the commutator of an operator A with H as ujA = [A, H]^ and the thermal average of A at 
temperature T as {A) = Tr[e~^/"^A]/Z with the partition function Z = Tr[e~^/"^]. The response function 
of "9 is introduced in the A:N x 47V matrix representation as K^ *t(i). Its {a^a') (g) (x,x') component is 
defined by 

K^%{t-x,x') ^ -t{[^Ut)A^^)x']h (4) 

Here, the Heisenberg operator for an operator A has been introduced as 

A{t) = e'"'Ae~'"K (5) 

The fermion anticommutation relations for Cx,s and cj,, ^, leads to 

^*,*t(0) = -i/ (6) 

with the AN x AN identity matrix I. The Fourier transform of eq. (|j) gives the susceptibility %^ ^t(w), 
which is nothing but the Green's function G{uj) in the 47V x 4A^ matrix representation. In more general, 
for arbitrary 4A^-component vector operators cj) and <p, their response function and susceptibilities can be 
defined in the 4A^ x 4A^ matrix representation as K^f^^{t) and X0(^t(w), respectively. Throughout this 



paper, the {a, a') (g) {x,x') component of the response function K^^^^{t) is defined by 

Kii,{t)^~z{\ct>%{t)A^^)i: 



(7) 



with the commutator and the anticommutator when 4> a-nd tp are operators with even and odd number of 
particles, respectively. The susceptibility X(t>,v\{'^) is given by the Fourier transform of i^^ ,^t(i); 



X0,yt(^) 



dte^'^'K^^^.it). 



(8) 



2.1 First- order profl^ction 
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As is well known JiiJ'13'' the first-order projection reproduces the Hartree-Fock theory. It also allows_iJS to 
construct the Dyson equation for the Green's function which is obtained from the perturbation theory Q* At 
this order of projection, only the zeroth- and first-order moments of the Green's functions appear explicitly. 
In the following, we give detailed arguments of the projection. 

The projection procedure is given by 

cj* = £(")* -I- to*, (9a) 

£(")^ii:^^^>,t(0)K-;^,(0), (9b) 

to* = (l-Pi)w*. (9c) 

Pi is defined as the first-order projection operator which operates to an arbitrary fermionic 4A'^-component 
vector operator (p as 



Pi0 = ii:^,^t(o)K^i^,(o)*. 

Pi satisfies the property of the projection operator, 

Pi=Pi. 
We note that the new operator Sui"^ satisfies 

-K'*,(5(i*)t(0) = -f^5(i*,*t(0) = 



(10) 

(11) 
(12) 



by definition. It allows us to express the single-particle Green's function G{ll}) = x^!t ift('^) in the form of a 
Dyson equation as 

G(cj) = i [g(°)-1(co) - i:i(co)l ' A:*,*t(o), 



The irreducible part of Xsui^ JSCd^'}'' i^) with respect to G^ '{cu) has been introduced; 

X5"*,(5i:j*)t(^) 

For the Hubbard model, one obtains 



(13) 

(14) 
(15) 



(16) 



^x2 = 4°L'P3CT0 - U6x^x' {Sx)pq(t - 



u 



(Ai 



V2 



■Pl<^a 



X 'Jx.x' 



(17a) 
(17b) 



Hereafter, the summation over i = 1 and 2 should be taken. We have introduced the local charge, spin and 



pairing operators at the atomic site x as 

nx = nx\ +nxi, (18a) 

Sx = cl-^a-ss'Cxs'/'^, (18b) 

^x = -j^{cx]Cxi + 4|4^), (18c) 

/\l = ^{cx^cxi-cl^cl^). (18d) 

The superconducting state with the isotropic s-wave pairing symmetry can be excluded in the repulsive 
Hubbard model (f/ > 0), i.e., (A^) = 0. The expectation value {Sx), if it does not vanish, is chosen to break 
the symmetry in the z-axis, i.e., {S^) — without loss of generality. Then, eq. (17a) and fc^ are reduced 
to 

£x,x' = 4°i'P3CT0 " f/(5^)5a;,a;'P0Cr3, (19a) 

SCj^x = U i -SnxP3Cro - SS^pQCrs j *a;, (19b) 

with Sux = Ux — (rix) and SS^ = 5| — (Sx)- 

If X'i(w) is ignored, then the present theory is reduced to the Hartree-Fock theory. The Hartree-Fock 
theory captures some of the many-body effects in the weak-coupling limit, while for the intermediate-to- 
strong couplings, the theoof breaks down. In fact, in the strong-coupling regime, the functional integral 
representation of the actiorO gives an expansion in a small parameter t/U around the Hartree-Fock solution. 
However, the Hartree-Fock theory fails to reproduce correct dynamical properties. In order to improve the 
Hartree-Fock theory, one must include the self-energy part corrections Si{u}). 

For this purpose, one can adopt perturbation expansions on a single-particle basis, one-loop approxima- 
tions, long-time or short-time approximations. However, these approximations do not give the energy- and 
momentum-dependent self-energy part in a correct manner. Particularly, if Si (uj) is self-consistenlhL,d£±ei:- 
mined by means of simple framework such as perturbation expansion or one- loop approxiinations,Cj'l£j'Ot3' 
the Mott insulator at half filling can not be reachedH3'E3 This is because arguments along this .line do not 
satisfy the Pauli principle and fail to give a correct high-energy behavior of the self-energy part .EH* 

The Pauli principle is a crucial aspect in reproducing a Mott insulator in the Hubbard model, since any 
deviation allows electrons (holes) to move through atomic sites already occupied by another electron (hole) 
with the same spin without any repulsion. Any incorrect treatment of the Pauli principle also violates the 
local moment sum ruleHI^ Only the zeroth- and first-order moments of the Green's function are correct 
in these theories, while the Pauli principle n^^ = Uxs enters a theory in the evaluation of the second- and 
higher-order moments. We note that the Pauli principle can be exactly satisfied in each order of the equations 
of motion. 

Another key factor to reproduce a Mott insulator is to correctly estimate the double occupancy. The 
double occupancy directly determines the local moment of the spin. In contrast with the self-consistent 
second-order perturbation theory, the dynamical mean-field theory supplemented with the iterative second- 
order perturbation theory reproduces the Mott insulator.H This is also because the Pauli principle and the 
double occupancy are correctly treated in solving the impurity problem. 

2.2 Second-order projection 

As is well known, the perturbation theory works more effectively for the self-energy part in the form of 
the Dyson equation rather than for the Green's function itself. Therefore, it is natural to consider the Dyson 
equation for the self-energy part. The operator projection theory is useful for this purpose. The second-order 
projection reproduces the two-pole approximation to the Green's function. It also allows us to construct the 
Dyson equation for the self-energy part. 

The second-order projection is given by 

wfc* ^ £(21)^ ^ £(22)^^2,^ + (5t2;(5w*, (20a) 

£(21) = K^s^^^^,{0)K^'^,{0), (20b) 



£(22) ^ fs:^,^^_(,^^)t (o)K-^^(,^^^t (0), (20c) 

Jwto* = (1 - F2)'i'5w*- (20d) 

P2 is defined as the second-order projection operator which operates to an arbitrary fermionic 4iV-coniponent 



vector operator as 



It satisfies the condition 



P2cP = K^MO)K-%,iO)^ 



K^js^j^ytiO'jK 



~i 



Sui'i/,{<^'^'^}^ 



(0)<5cj*. 



P9 -P2. 



From eq. (|2C|), one obtains the irreducible self-energy part 






n(0)-l 



to 



)-i:2M 



=-(21) 



CjJ — £ 
irr 



(22) 



(21) 

(22) 

(23a) 

(23b) 
(23c) 



-S'2(w) — -''X!5ti5i.*,(<5'i'5'i*)t('^)-^<5^*,(5LJ*)t('^) 

Here, we have introduced the irreducible part oi XsiiiSiii^JSdjSdj^yti^) with respect to G^ '{cu) and ^1(0;) as 

?^(5ii(5ii*,(<5ii'5'^*)t '^-' 

For the repulsive Hubbard model, 

'1 



(24) 



^^a: = (( -^^nxpscro - SS^po<73 



■x^x' = -tx%P-i^o - (M2P3CT0 - {SDpqct^) 5, 



-x,x 
.(22) 



'x.x' ^ 



(25a) 
(25b) 
(25c) 



*af,x' = tx,x' {SCxPsSCx' P3)M, 

Oi-a; ^ •'-a; (•'-a; 
1 



-1 



1 



Cx = -SuxPscra - SSx ■ Poo- + —^Axpjao, 

M2 = P2 P0<^0 + [(1 - ("-))ekin + £cor2] Mx^ , 



(0) 
Ai2 



: p-U (l- -(ria 



£kin = ^77 X!*fc ( 2 ^"'0;'^^^3^" ~ ('5'o;fe)PoO-3 



iV 



with 



s 
^<l,k = 2 / .%,-q/?...^'^ss'Ck+k/2,s'- 



(25d) 
(25e) 

(25f) 

(25g) 
(25h) 

(25i) 

(26a) 
(26b) 



The kinetic energy per site summed over spins is expressed as i^kin — ^Trekin(po+P3) — —jj X]fc s ^ki'^ks'^ks) ■ 
Below, a momentum independent energy shift ecor2 due to two-site correlated hopping terms which take the 
form tx.x{cxsCxs'Snx-s) in eq. (25c) are neglected. We note that in the particle-hole symmetric case, ecor2 



vanishes. Then, in the case of (S^) — 0, the remaining operator becomes 

-'Utx,x6CxP3'i'x- (27) 

A general two-pole approximation by neglecting E2o('-^) and its approximate treatment will be given in 



2.3 and § 2.4, respectively. 



2.3 Neglecting S2{^) — Two-pole approximation 

If X'2(w) and-jbhe momentum-dependence of £(^2) g^j-g neglected, the present formalism gives the Hubbard I 
approximation^) The momentum dependence of e^^^' appears through the equal-time short-ranged charge, 
spin and isotropic s-wave pairing correlations. XWs modifies the dispersions of the separated Hubbard bands. 
This yields the Roth's two-pole approximation^ 

The two-pole approximation has been found to give a good starting point in discussing main single-particle 
dispersions near the MIT. Two dispersions of the energetically sei)|ai;ated Hubbard bands obtained with this 
approximation resemble those obtained by QMC calculations, t^lL^P while it does not reproduca. tiLe-^elf- 

are also known to violate the Luttinger sum rule due to the formation 
of two distinct Hubbard bands. To overcome these problems, it is required to include the higher-order 
dynamics. 

Here, we briefly discuss a general two-pole approximation. In the symmetry-unbroken phase, the electron 
Green's function with momentum k and spin s reads 



energy part at cj = obtained byJCanamori in the low density limit. Two-pole approximationstallaEj-' as 
well as Hubbard approximationsB^ E ' 



Gc,s(w,fc) 



(0) (|)^(n)(2^(n)) 



(28) 



Here, the (1,1) components of eqs. (25c) and (25d) and ecor2 have been introduced in the momentum space 
as 

(0) 2(l-(n))gMn+4gcor2 ._. 

'^ = -'^-^'^ (n)(2-(n)) ' (29) 

tf^ and Ecoi-2, respectively. The above Green's function has two poles at w = w^^ with their residues z, 
where 

'-fe=(4"^+^~fe±Ar)/2, (30) 

.(0) ~. \/aMH 



z^^il±{el"-e„)/A^r)/2. (31) 



Here, we have introduced 



A 



r = v'c/2(n>(2-(n)) + (4°)-efc)2, (32) 

which generate the direct energy gap at the momentum k between the two poles lying in the upper and lower 
Hubbard bands. The momentum dependence diminishes as U/t diverges. It is easily found that when the 
collective degrees of freedom are completely frozen, i.e., i^ = 0, and the third term in eq. (E9h is neglected, 
then eq. (Eq) recovers the Hubbard I solution. 

Next, we consider the case where equal-time two-particle correlations with the momentum Q are significant. 
Then, through the growth of short-range correlations in eq. (25d), ik behaves like tk+Q so that eq. ( |2^ ) 
continuously connects two different phases that preserves and breaks the translational symmetry. Actually, 
when the commensurate AF long-ranged order is stabilized at the ground state in the simple square-lattice 
Hubbard model, the nonzero ordered moment in this case creates new modes in single-particle excitation 
spectra. The new modes emerge due to the spontaneous breaking of the translational symmetry, namely 
the folding of the Brillouin zone. Therefore, when the higher-order dynamics is ignored, a discontinuous AF 
transition can occur, as it does in the Hartree-Fock approximation. However, this is not always valid. Upon 
continuous phase transition, this discontinuous feature should be removed by a gradual growth of the AF 
shadows. Any new mode leading to a continuous phase transition does not appear discontinuously at the 
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phase transition, if and only if the condition given by i;^ = tfc+o or cquivalcntly 






for the different sublattices 
for the same sublattice 



is satisfied at T = and tx^x' changes as a continuous function of temperature. So far, few theoriesJiaMe 
reproduced such shadow-band effect even quahtatively within the one- loop approximations to X'i(a')£3'Ei|) 
The two-pole approximation can capture the shadows and continuously approach the AF phase from the 
paramagnetic phase. In this case, however, the two-pole approximation only gives the SDW bands similar 
to the Hartree-Fock results. Apart from the shadow features, the approximation does not introduce any 
additional feature beyond the Hartree-Fock approximation. In § 0, it will turn out that low-energy narrow 
bands playing a role of quasiparticles do emerge as significant additional features. 

We note that as is clear from eq. (25d), ifc includes the superexchange term. This suggests that the 
momentum dependence of the self-energy part appears through the superexchange term. For further im- 
provements of the theory, higher-order dynamics contained in £'2(0;) has to be considered. The following 
arguments on the expectation value of the Hamiltonian illustrate this necessity. 

In general, the expectation value of the Hamiltonian is given by 

U{n) 



(H)=T 



'n,fe-S 



-ik 



■'lc,s 



(w«,fc) 



Gc^s(^n,fc) 



(33) 



In the two-pole approximation without I]2{^)^ we obtain 



m 



k,± 



u 



± 



{n){2-{n)) 



(34) 



as a function of ti^. Here, we have introduced a fermionic Matsubara frequency a;„ = (2n -I- 1)t:T and 



the Fermi distribution function f{uj) — [e^l^ 
component of J7i(w). At (n) = 1, we obtain 



1) . Sie,s(cL;, fc) denotes the momentum-fc spin-s electron 



(^> 



fc 



[4 



(ife+^fc)74]/C/ + 0(iVC^') 



(35) 



when U '^ T, t. This clarifies that the superexchange term appears through eq. (25d). However, there exist 
extra terms of order of t'^ /U in eq. (pq) in addition to the superexchange term. Since the term of order of 
t'^ /U should be uniquely determined in terms of the strong-coupling expansion, the additional terms are to 
be cancelled with the terms from S2{^)- 



2.4 Higher- order dynamics 

Here, we consider the second-order self-energy part 1^2 (^) defined by eq. (23c). In principle, one can 
continue the projection procedure to derive the Dyson equation for 172 (tj). The results of the third-order 
projection will be published elsewhere. In this paper, we employ an approximation to obtain ^72(0;). In the 
symmetry-unbroken phase, its off-diagonal components vanish and it does not depend on the spin index. 
Then, the spin indices can be dropped out in the following notations and we have only to discuss the 
spin-independent electron component expressed as 

^^'=^'^' ^^ " (n)(2-(n)) ^^'"^""^'('^"'^"'^)' ^'^' ^^' ^^^^ 

Several approximate methods are available to calculate I]2o(^, fc): the perturbation theory, one- loop ap- 
proximation with two-particle susceptibilities obtained by other methods such as the_aelf-consistent RPA,^^ 
the two-particle self-consistent methodE^'Ed-' or the extension of the Roth's methodH by Mancini et alE3> 
which exploits equations of motion of the response functions of "^x with respect to Uxs'^x' ■ Here, at least, 
the following two points are required in calculating E2c(^, k) along the present approach: One is that they 
should give a correct estimate of the moment of Ti2c{'-^,k), thus the fourth-order moment of the Green's 
function. Th£ ijtlicr_is_that the short-ranged correlations which are neglected in previous equation-of- motion 
approachesBDQ'Ej'EErH' should be taken into consideration. 

The perturbation theory gives an approximately good estimate of the moment of S2c('.i^, k). However, it 
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ignores the short-ranged correlations. Since one of our main aims is to take account of both the short-ranged 
correlations and the Mott-Hubbard picture due to the strong local repulsion, the perturbation theory is 
not sufficient. A possible prescription to take account of the both is to decouple the three-body correlation 
function S2e(w,fc) into the convolution of the single-particle Green's function and two-particle dynamical 
susceptibilities. l— . 

We note that ^2(^) has been evaluated within a two-site level by means of a different decoupling scheme.EJ 
However, short-ranged correlations are ignored there. If the AF short-ranged spin correlations are taken into 
account, the AF shadow bands of the Hubbard bands should appear, as we will obtain below. 

First, we give our decoupling scheme. In the symmetry-unbroken phase, the spin-independent electron 
component of the second-order electronic self-energy part is evaluated by means of the following decoupling 
approximation as 

T.2c{(^n,k) = j-^j^ j-^ ^ T.2co{'^n,k), (37a) 



(n>(2-(n)) 



0=n,SA 



S2oo(w„,fc) = — ^Ao(A;;q)^Go(cj„ -f]™,^-^) 

m,q 

X Xo,ot(ri„,q) 
for O = n and 5', 

--^^AA(fc;9)^Go(^m-w„, 

m,q 

X XA,At(f^m,g), 



S2eA(w„,fc) 



k) 



(37b) 



(37c) 



where flm = 2T:mT is a bosonic Matsubara frequency, Ge is the spin-independent electron part of the 
Green's function and Xn,n, Xs\S' s-nd XA.At are the charge, spin and local-pair susceptibilities, respectively. 
The following vertices have also been introduced; 



with 



A„(fc;q) = -Tfc.q, 
As4fc;q) = -2Tfc.q, 

Asx.„(fc;qr) = -2ifc_q, 

AA(fc;9) = -2ife_q, 
2(1 -(n)) 



^k.q ''k—q ^k 



(n)(2-(n)) 



£^kin- 



(38a) 
(38b) 
(38c) 
(38d) 

(39) 



Here, tf^ is the (1, 1) component oi t\ . These equations are straightforwardly derived from eqs. (23c) and 
(27). 

Figure 1 shows the Feynman diagrams of the present decoupling approximation to the above five (n, 5^, 
5^, S^ and A) components of £20 in the imaginary time and the real space. Here, the vertices have the 
following representations in the real space; 



A„(a;;a;i,a;2) = Sxi,x 

As-{x;xi,X2) = 26xuX2 
-26.. 

As^-y{x;Xi,X2) 



■ 2(1 - (ri))i?kin ^ ,(22) ' 

'^"■"^ {n){2~{n)) + "^"\ 



■ 2(l-(n))gkin ^ ,(22) 

^'=^^ (n)(2-(n)) + =^'=^\ 



,t 



X,Xil'Xi,X2 ) 



2dx.XiL'Xi,X2 5 



AA(a;;a;i,a;2) = -25x^xitxj^,x2- 



(40a) 

(40b) 
(40c) 
(40d) 



12 



X„,„(t'-t,x,'-x,) 






5:2e.(T-T';X,X') = A'^'"" ^ '''"''A 



(a) 



An(x';xi',X2') An(x;xi,X2) 

Zs'.s'Ct'-cxi'-xi) 
.. < 



S2eS'(T-x';X,X') =^ 

(b) 



■ ■^ _ ; (i',x,') (t,x,) ^ 



As'(x';xi',X2') As'(x;xi,X2) 

ZA-A*(T:-t',Xi-Xi') 



/ \ 



(C) 



Aa(x';xi',X2') A4(x;xi,X2) 



Fig. 1. Diagrammatic representations of the second-order self-energy part in our decoupling scheme. It can be decomposed 
into three parts shown in (a), (b) and (c) divided by {n){2 — (n)), as in eq. (37a). In (b), the superscript i takes x, y and z. 
Details are explained in the text. 



From eq. (^7|), it turns out that the second-order self-energy part takes the similar expression to the 
usual self-energy part, if one takes decoupling approximations to evaluate them. The difference between 
the self-energy part and the second-order self-energy part appears in their coupling constants. The coupling 
constant for a decoupled form of the self-energy part is C/, while that for the second-order self-energy part 
is Ae)(fe; q), of order of t. 

The momentum dependences of h-o{k] q) determine detailed contributions of two-particle correlations to 
E2e(w, fc). To discuss these contributions, for simplicity, we consider the Hubbard model with the electron 
transfers restricted to the nearest neighbors. When AF spin correlations are so strong that the equal-time 
nearest-neighbor spin correlation has a large negative value, it happens that i^ approximates ifc+n = — ifc 
with the staggered wave vector Q. In such case, the first and second terms in eq. (|3^) nearly cancel 
each other for momentum q ~ Q. Particularly, at half filling, eq. (^) suggests that T^.g nearly vanishes. 
Therefore, when AF spin correlations develop at half filling, the main contributions to S2e((j-', A;) are those 
from E2e5=:,B(w, fc) and E2eA('^,fc)- In the strong-coupling case, the double occupancy is strongly reduced. 
This suppresses the contribution from E2eA(w,A;). Finally, eqs. (37b), ( ]3^ ) and ( p9| ) suggest that for the 
momenta lying on the magnetic Brillouin zone boundary where tj^ and thus tfc_n vanish, all the contributions 
from the two-particle fluctuations with their momenta q = Q and Q to S]2c(w,fc) disappear at half filling. 
Even away from half filling, only the charge and the longitudinal fluctuations dominantly contribute to 
I]2c(w, k). On the other hand, away from the magnetic Brillouin zone boundary, E2o(^, k) is dominated by 
the contributions from the transverse AF spin fluctuations. Therefore, modiflcations of dispersions from the 
Hubbard bands stem from the spin excitation properties under strong AF correlations. In § |^, we actually 
see additional low-energy bands develop near half filling at low temperatures. We will discuss the issue in 

§1 

Next, we must obtain the two-particle susceptibilities in an appropriate method. As we have already 
emphasized, effects of the short-ranged correlations must be incorporated into the present theory for a better 
description of the strongly correlated electrons. For this purpose, one should include the Stoner enhancement 
factor resulting from a summation of ladder or bubble diagrams and also obtained from the Gaussian 
approximation in the functional integral techniqueHj^ Although one can easily calculate the dynamical 
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spin susceptibilities using the method developed by Mancini et aL,c3> such enhancement factor is not taken 
into account in the approximation. Another requirement already noted above is that in the framework 
of the present 0PM, the-higher-order dynamics accounted by S2{'^) raust reproduce its correct moment. 
The self-consistent RPAEJ'EJ) is based on the equation of motion of two-particle operators. Therefore, it 
seems that the approximation may be appropriate. The self-consistent RPA, however, does not satisfy the 
local moment sum rules explicitly given below. Instead of these approximate methods, the two-particle 
self-consistent (TPSC) method determines the irreducible vertices in the RPA-like forms of the two-particle 
susceptibilities so that the lacal-moment sum rules are fulfilled for the ±wo-particle properties, not only 
charge and spin correlational^ but also local s-wave pairing correlations.EJ In this paper, we employ the 
TPSC susceptibilities as in the previous paper. If one employs the perturbation theory or the self-consistent 
RPA, they also reproduce a Mott insulator. However, they give weaker AF spin correlations than the TPSC 
method. 

Within the two-particle self-consistent approximationj£3' EJ-* the susceptibilities are calculated from 

XS^M^rn,k)= '^^^^"7„^^.. , (41b) 

1 - i sXph("m,«;) 

XA,At(^™,fc) - /^Pp("-fe) (41c) 

with an approximation Fs = 'HJ {nx^Uxi) / {n) and the self-consistency conditions for Fg, Fc and Fp given 
by the local moment sum rules; 



T 



'^^ Xn,n{^m,k) = {n) +2{nx]nxi) - {n) , (42a) 



N ^^ 

m.k 

j;^'^Xs\s4^n,,k) = (n) - 2(na:^na:i), (42b) 

m.k 

-^ X! ^^-^^ (^"' ^^ " {"^x^rixi) ■ (42c) 

ni,k 

Xph and Xpp are the particle-hole and particle-particle susceptibilities calculated from the bare Green's 
functions. 

§3. Numerical results 

In the previous letter ,El it was shown that the filling-control MIT is well described by the present method. 
In this section, we calculate the single-particle density of states, spectral functions, dispersions, momentum 
distributions and the chemical potential shift. Their dependences on T, {n) and/or U are also discussed in 
detail. 

The following results have been obtained with 2048 Matsubara frequencies in the 32 x 32 lattice for two- 
dimensional systems and in the 128_aites chain for one-dimensional systems. For analytic continuations, we 
have used the Fade approximation.ElP 

3.1 Mott-insulating properties at half filling in two dimensions 

First, we give results at half filling. Results on the local density of states p{lo) = —{ttN)~^ ^^ ImGe(w, A;) 
for U/t = 1, 4 and 8 are shown in Fig. 2 (a), (b) and (c), respectively, at several temperatures. At 
temperatures low enough, a Mott-insulating gap grows. This indicates that the ground states are Mott 
insulators. At a much smaller coupling U/t = 0.2, we have also obtained a Mott gap. Then, we expect 
that the present method reproduces a Mott insulator at any finite value of U. This Mott-insulating feature 
is accompanied by sharp ftK^ks in p{uj) at the low-energy edges of the Hubbard bands in contrast with the 
Hubbard approximations. tlB Excitation spectra around (tt, 0) and (Or-fr) mainly contribute to these sharp 
peaks. These sharp peaks have also been found in QMC calculationsrS-' As the temperature increases, the 
spectral weights fill in the gap and the intensities of the peaks diminish. At high temperatures T > U/2, 
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a low-energy incoherent band appears inside the Mott-Hubbard pseudogap separating the Hubbard bands. 
This is particularly prominent in the strong-coupling regime. The whole band structure becomes similar to 
the free-electron tight-binding band, particularly in the weak-coupling regime. The excitations in this band 
are incoherent and can not be described by well-defined quasiparticles. This is confirmed by the smooth 
momentum distributions n^ in the momentum space, as shown in Fig. 3 (c). 

In the Fermi liquids, in general, momentum distributions have jumps along the Fermi surface. On the 
other hand, in the present case, i.e., the repulsive Hubbard model on the square lattice at half filling, the 
momentum distribution exhibits neither a jump nor a singularity. It is a smooth function of the momentum, 
as shown in Fig. 3. This means that the effective Fermi temperature vanishes in this case. This is naturally 
expected at low temperatures, since the system is a Mott insulator at T = 0. At high temperatures, low- 
energy excitations are recovered around lo — 0. In spite of such metallic features, the momentum dependence 
of the momentum distributions monotonically decreases. This means that the system approaches the classical 
limit where n/, is independent of fc. Similar situations with a suppressed Fermi temperature also occur in 



the doped cases, as discussed in § ^^. Figure 3 also shows that the momentum dependence of n^ decreases 
with increasing U. The deviations of n^ from the line for U — become smaller for higher temperatures. 
This is naturally expected, since these states at T > are adiabatically connected with the free-electron 
system. 

To clarify an accuracy of the present method, the T — QMC results for U/t — 4 obtained by Assaad and 
ImadaLiP are also shown for comparison in Fig. 3 (b). Those for U/t = 8 by AssaacO* are shown in Fig. 3 
(c). Compared with QMC calculations, the error turns out to be less than 10% for U/t = 4. The error is 
even smaller for U/t — 8 than for U/t — 4. This small discrepancy originates from an overestimate of the 
Mott gap for smaller values of U. Though the Mott gap is actually reduced from U due to fluctuations, 
the reduction ratio is still small compared with QMC calculations. The discrepancy in the evaluation of 
the Mott gap will be examined later by showing the single-particle dispersions as well as by obtaining the 



chemical potential as a function of the doping concentration in § 3.2. 

The energy and momentum dependences of these single-particle excitations can be understood by calcu- 
lating the spectral function, 

A{oj,k) = --ImGc{uJ,k). (43) 

n 

Figure 4 shows the single-particle spectral function A{uj, k) for U/t = 4 at T/t = 0.02 (a), 0.5 (b) and 4 (c). 
The upper and lower Hubbard bands nearly separated by the Mott-Hubbard pseudogap exist at T/t = 0.5 as 
shown in Fig. 4 (b). At lower temperatures, the growth of AF spin fluctuations produces their shadow bands 
due to the folding of the Brillouin zone in the AF ground state. This results in a four-band structure, as 
shown in Fig. 4 (a) for T/t = 0.02. This evolution of the single-particle dispersions from the Hubbard bands 
into such four-baiLd_structure at low temperatures T < J has also been obtained with QMC calculations 
by Grober et a/EJ'LJ We will discuss the dispersions in more detail later. On the other hand, at such 
higher temperatures as T > U/2, the Mott-Hubbard pseudogap disappears and the low-energy incoherent 
band develops, as shown in Fig. 4 (c). It seems that this realizes a similar band structure to a free-electron 
tight-binding dispersion. However, as we mentioned above, this band does not introduce any abrupt change 
in the momentum distribution. This reflects that the Fermi temperature vanishes for this half-filled case. 
Therefore, for arbitrary values of C/ > and T, it is not justified to describe the electronic structures 
of the simple two-dimensional half-filled Hubbard model in terms of the Fermi-liquid theory. Though these 
incoherent metallic states at high temperatures are adiabatically connected with a free-electron ground state, 
it does not necessarily mean that a simple weak-coupling perturbation theory effectively works for describing 
Mott-insulating properties at low temperatures. Only those theories which describe both the Fermi liquid 
and the non-Fermi liquid reproduce correct electronic properties. 

Lastly, we discuss an evolution of the electron dispersions at U/t = 4 as a function of temperature. 
Figure 5 shows the peak positions in A{uj, k) at the same temperatures as given in Fig. 4. The QMC 
results at T = (& are also given as the circles in Fig. 5 (a), while the present results are marked as the 
squares. This figure indicates the overall similarity of the present results to the QMC results, except for an 
overestimate of the Mott gap. It should be noted that the dispersion around (7r,0) and (0,7r) obtained in 
the presftnt -calculations is much ffatter than the IfabtjacdEP I approximation, in agreement with the QMC 
resultsEaE^I'Eil'ESP and the two-pole approximationJla'td'H Moreover, the four-band structure composed of 
the upper and lower Hubbard bands and their AF shadows takes the form of a superposition of tato-SDW-like 
bands and two low-energy narrow bands. Such structure is consistent with the QMC results.E^EJ As the 
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temperature increases above J, strong short-ranged AF spin correlations disappear. This smears out the AF 
shadows of the Hubbard bands. Then, the main structure is characterized by the Hubbard bands, as shown 
in Fig. 4 (b). At higher temperatures T > U/2, an extra low-energy narrow band crossing the chemical 
potential appears without any coherence, as already mentioned and also shown in Fig. 5 (c). The main 
structure shown in Fig. 4 (c) resembles the free-electron tight-binding band. 

3.2 Metallic properties in doped cases 

Next, we give the results for doped cases with (n) — 0.87. Since electron transfers are restricted to the 
nearest neighbors, the following results also hold for (n) = 1.13 after particle-hole transformation. 

Results of the local density of states p{uj) in the case of (n) = 0.87 are shown in IriS|-§ for U/t = 1 (a), 4 (b) 
and 8 (c) at several choices of temperatures. As in the two-pole approximations,t3E3' the upper and lower 
Hubbard bands are nearly separated. In addition, weak incoherent spectra remain between the Hubbaixl 
bands even at T/t — 0.02. For U/t — 8, these features have also been obtained in QMC calculations. H 
For all the used values of U, jfl-sharp peak develops near the Fermi level as the temperature decreases, in 
agreement with QMC results.ES This sharp peak emerging at low enough tempemtures is also reminiscent 
of the paramagnetic solution of the doped Hubbard model in infinite dimensions .ES) However, an important 
difference here from the infinite dimensions is that the sharp peak is pinned at the Hubbard gap edge. At 
high temperatures T > U/2, A{u!, k) exhibits a rather broad peak between the Hubbard bands, as it does 
at half filling. 

Figure 7 shows the momentum distribution n^ for the same electron filling as in Fig. 6, (n) ~ 0.87 and 
U/t = 1 (a), 4 (b) and 8 (c) at various temperatures. In the weak-coupling regime, the results for U/t — 1 
and T/t = 0.02 suggest that the Fermi surface at which the momentum distribution jumps at T = exists 
and that it appears to agree with the free-electron case, as shown in Fig. 7 (a). This is consistent with the 
statement that this system belongs to a Fermi liquid. The difference in n^. between U/t = 1 and cases lies 
in the value of the Fermi jump. On the other hand, for U/t = 4 and 8, Fig. 7 (b) and (c) show that the 
momentum distribution has the maximum slope outside the magnetic Brillouin zone along both directions 
from (tt, 0) to (7r,7r) and from (0,0) to (tTjTt). This suggests that the obtained "Fermi surface" is larger 
than the Luttinger volume at least at T/t = 0.02. The "Fermi jump" determined at T/t = 0.02 decreases 
with increasing U. This leads to the smaller momentum dependence for larger values of U. At higher 
temperatures, the "Fermi jumps" are smeared out. At T/t < U/2, the slopes of n^ at the "Fermi surface" 
appear to be smaller than in the free-electron cases. This indicates that the characteristic temperature 
for smearing the "jumps" in the momentum distribution at the "Fermi surface" decreases with increase in 
U. This occurs because the coherence in single-particle excitations is suppressed due to correlation effects 
near the MIT. At higher temperatures, n^ approximates the free-electron value. With further increasing 
temperature, the quantum-mechanical degeneracy of electrons disappears. Then, the momentum dependence 
of n/j gradually disappears. 

There is another remarkable property in the momentum distribution n^. Figures 7 (a) and (b) show that 
at T/t — 0.02, after the sudden decrease in nj~ at the "Fermi surface", it exhibits an increase along both 
(0, 0)-(7r, tt) and (0, 0)-(7r, 0)-(7r, tt). If this second large change in Uf, evolves into another jump at T = 0, 
this violates the Luttinger sum rule. 

One main reason that we have obtained the Fermi volume inconsistent with the Luttinger volume at 
T/t — 0.02 is the downward shift of the main dispersion in the lower Hubbard band around the momentum 
(7r,7r). As shown for U/t = 8 in Figs. 8 and 9 and discussed later, the low-energy single-particle dispersion 
in the lower Hubbard band is significantly weak outside the magnetic Brillouin zone boundary as well as 
around the (tt, 0) and (0, tt) momenta. Figure 8 also shows that spectral weights in the lower Hubbard band 
outside the magnetic Brillouin zone boundary are strongly reduced. This makes it difficult to speculate 
the T = properties from the present results. More detailed and rigorous considerations on the Luttinger 
sum rule requires results at T = 0. In addition, the dispersion and the momentum distribution are still 
sensitive to a choice of approximations: Matsumoto and Mancini obtained a large Fermi vnlume for similar 
choices of parameters by employing a different decoupling scheme based on a local picturcEHF It is not clear 
whether their method satisfies the Luttinger sum rule as the system approaches the half filling. Particularly, 
detailed low-energy structures depend on the spin fluctuation properties. Further studies are required for 
this problem. 

The single-particle spectral functions and the corresponding dispersions for the filling (n) = 0.87 and 
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U /t = 8 are given in Figs. 8 and 9. The energetically separated Hubbard bands clearly appear even in this 
doped case at all the values of temperature we have used. The main structures are similar to the QMC 
results.Ea^ Particularly, the single-particle dispersion in the lower Hubbard band becomes strongly weak not 
only outside the magnetic Brjllpuin zone boundary but also around the (7r_Q\_and (0, tt) momenta. This 
agrees with QMC calculation£3^ and results of the two-pole approximation.tSEiJ At T/t — 0.02, an extra 
low-energy band appears at the top of the lower Hubbard band and is pinned at the top. This seems to evolve 
into the low-energy band obtained as an AF shadow of the lower Hubbard band at half filling. Another 
important feature is, as mentioned above, that the deviation from the Luttinger sum rule is found at least 
for the temperatures investigated here. Actually, QMC calculations show that for T/t ~ 0.5, the (7r,0) level 
is below the Fermi level.H' In Fig. 10, we also show the results of A{lj^ k) for the same filling (n) = 0.87 but 
with U /t = 4 at T/t — 0.02, 0.5 and 4. They appears to be similar to those with U /t = 8, except a smaller 
value of the Hubbard gap. 

Next, we consider the chemical potential compared with the particle-hole symmetric position /i — C//2 versus 
the electron filling (n) shown in Fig. 11. Near half filling, thn^/i-(n) curve shows that k — d{n)/dfj, appears 
to diverge as (n) ^- 1 in agreement with the QMC results.H For U/t — 4 and 8, the chemical potential 
shows strong (n) dependence for (n) very close to 1, namely for (ji) > 0.95 and 0.9, respectively. This occurs 
because slight but finite spectral weights still remain inside the zero-temperature Mott- Hubbard gap of order 
of U due to finite-temperature effects. Because of the smaller residual states for larger U, fx dependence 
of n becomes larger for larger values of U . In addition, with increase in U , the difference in momentum 
distributions between the momenta (0, 0) and (tt, tt) decreases, due to larger low-energy spectral weights 
around (tt, tt) and (0, 0) in the lower and the upper Hubbard bands, respectively. They enlarge the range 
of the hole concentrations where the finite-temperature effect appears. At strictly T = 0, d{n) /dji — > oo is 
expected with {n) — > 1. We note that the divergence of the compressibility is in contrast with the prediction 
of the Brinkman-Rice picture where the compressibility remains finite. This important difference is more 
clearly seen in the density of states at half filling in Fig. 2. The nearly diverging density of states near the edge 
of the lower Hubbard band appears due to the low-energy band with extremely weak dispersion especially 
around (tt, 0) and (0,7r). Such low-energy band near the level of (7r,0) and (0,7r) momenta is beyond the 
Brinkman-Rice picture. As we will argue in § 0, strong AF spin correlations and their contributions to the 
strong momentum dependence of the self-energy part produce such sharp contrast with the Brinkman-Rice 
picture and the dynamical mean-field theory. 

Lastly, we discuss the Fermi velocity and damping rate that characterizes the low-energy single-particle 
exciations. Figure 12 shows the Fermi velocities up = [|Ve*(fc)|]fe=fcp and the damping rates 7^^ at the Fermi 
momenta along both (0, 0)-(7r,0)-(7r, tt) and (0, 0)-(7r,7r). It shows the divergence of the quasiparticle mass 
in both momentum regions near (7r,0) and (7r/2,7r/2). although the anisotropy is large. The diverging mass 
itself is consistent with the Brinkman-Rice picturcEj-* These results reflect the generic suppression of the 
Fermi degeneracy near the MIT. However, the enhanced density of states due to the low-energy band is pinned 
at the top of the lower Hubbard band. This is the origin of new physics not contained in the Brinkman-Rice 
picture. Furthermore, quasiparticle damping rates given as |ImEie(w = 0, fcp)! exhibit a prominent increase 
as (n) -^ 0. Thefore, quasiparticles become ill-defined at least at T = 0.02t, as one approaches the MIT. In 
fact, for 0.94 < (n) < 1, low-energy excitations near the Fermi level and momenta in Aiuj^k) exhibit a too 
broad structure to recognize quasiparticles with well-defined velocities and damping rates. Therefore, the 
Fermi surface is also ill-defined in this filling region at least at this temperature T/t = 0.02. Further studies 
on the low-energy single-particle excitations are required in order to judge the validity of the Fermi-liquid 
theory. 



3.3 Mott insulator in one dimension 

As mentioned in § y, one-electron £xcitation spectra of the Hubbard model in one dimension have been 
understood in detail only for U -^ 00 J3 where spinon and holons are decoupled at an arbitrary energy scale. 
This produces spinon and holon branches in one-electron excitations. These branches do not have poles due 
to the one-dimensionality, leading to the Luttinger liquid instead of the Fermi liquid. In the case oi U —^ 00 
at half filling, it is only at |fc| < 7r/2 that the spinon branch has finite spectral weights below the chemical 
potential because of the conservation laws for energy and momentum in creating a holon and annihilating 
a spinon. For the half-filled Hubbard model with a finite value of [/, their one-electron spectral properties 
are still not clear, since the spinous and the holons can not be decoupled at all energies. In this case, QMC 
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calculationala^ give a similarity to the SBMF predictioncj) and predict that at low but finite temperatures 
without any AF long-ranged order, the shadow structures develop. Here, we apply the present formalism to 
the one-dimensional Hubbard model and compare the results in two dimensions with those in one dimension. 

Figure 13 shows our results for the spectral function A{uj, k) of the one-dimensional Hubbard model with 
U /t = 4 at half filling. Similarly to the two-dimensional case, a Mott gap opens and separates the Hubbard 
bands. At T/t = 0.2, single-particle excitations at the momenta k = ±7r/2 become the closest to the Fermi 
level. This refiects that a negative value of t/t provided by AF short-ranged spin correlations yields SDW- 
like bands instead of the simple Hubbard bands, in terms of the two-pole approxiniation. Apart from the 
overestimate of the Mott gap, these main structures are similar to the QMC results.© 

At lower temperature T/t = 0.01, the four-band structure emerges in the spectral function, as in two 
dimensions. The origin of the four-band structure is the same as that in the two-dimensional case: The 
AF shadows of the Hubbard bands develop due to the growth of AF short-ranged spin fiuctuations. The 
structure composed of the Hubbard bands and their shadows yields the four bands. In the low-energy band 
emerging at the top (bottom) of the lower (upper) Hubbard band, spectral weights appear not onJx for 
|fc| < 7r/2 (|fc| > 7r/2) but also for |fc| > 7r/2 (i^r^ 7r/2). This contrasts with the result a,t U ^ cx)B or 
exact diagonalization results in the t-J model JI3' 13 which show that the spinon branch appearing at the 
top of the lower Hubbard band has finite spectral weights only for |fc| < 7r/2. The present results for an 
intermediate coupling can not be well understood from the arguments at L/ — > oo in terms of the decoupled 
spinous and holons. This indicates strong mixing effects between spinous and holons. However, these low- 
energy bands may still be sensitive to the choice of an approximation to £20 (w, k), as previously mentioned. 
Therefore, further studies are required to clarify these low-energy bands in more detail. 

§4. Total View of the Mott Transition in Tvi^o Dimensions 

In this section, we discuss the route to the Mott metal-insulator transition (MIT) in the two-dimensional 
Hubbard model. We provide a unified view of the Mott-Hubbard, Brinkman-Rice and Slater pictures. Then, 
we clarify to what extent these old pictures hold and how they are modified. 

First, we note that as one approaches the MIT, the local spin moments develops as expressed in eq. (42b). 
It increases with increasing U through the reduction of the double occupanpiv:.. |This gives rise to the splitting 
of the upper and lower Hubbard bands as in the Mott-Hubbard picturcBu'EP At low temperatures, this 
large local spin moments tend to align antiferromagnetically to gain a superexchange energy, as in the Slater 
picture. til' Then, these local AF spin correlations contribute to formation of the SDW-like dispersions in the 
Hubbard bands through the momentum dependence of the self-energy part as |ifl-|the two-pole approximation 
at half filling.E3 We note that not only thejSimple Hubbard approximationsoQ-EP but also the paramagnetic 
solution of the dynamical mean-field theoryH gives a simple cosine band in each Hubbard band but not the 
SDW-like dispersion. This is because the momentum dependence of the self-energy part is oversimplified 
or neglected. At low temperatures, a narrow hand appears at the top of the lower Hubbard band for hole- 
doped cases, as in the Brinkman-Rice pictured However, the origin of the narrow band at the gap edge is 
clearly different from the Brinkman-Rice and dynamical mean-field theories, because the enhanced density 
of states already exists at the gap edge even in the insulator phase. Since the dispersions in the Hubbard 
bands take the form of SDW-like dispersions, the level of the narrow band is naturally located near the 
magnetic Brillouin zone but not (tt, tt) in contrast with the Brinkman-Rice picture. The chemical potential 
shifts within the narrow band at low doping concentrations. The velocity of these low-energy single-particle 
excitations increases with increasing the doping concentration. But the level of the narrow band compared 
with the SDW-like dispersion in the lower Hubbard band does not change and pinned at the edge of the 
Hubbard band, in contrast with the results of the dynamical mean-field theory. These lead to the divergence 
of the charge compressibility. 

Further studies on the low-energy narrow band are necessary not only to discuss the Luttinger sum rule 
but also to clarify the mechanism of the high-Tc cuprate superconductivity, since it seems to play a significant 
role in recovering the quantum-mechanical coherence. 

Next, we summarize the modifications of the three old pictures of the MIT clarified by the present theory. 

The Mott-Hubbard picture holds at high temperatures where the spin correlations are random. However, 
when the spin correlations develop at low temperatures, this picture breaks down. Instead, equal-time local 
AF spin correlations yield the SDW-like dispersions in the Hubbard bands, as in the two-pole approximation 
at half filling. Simultaneously, the original disperion of the lower (upper) Hubbard band shrinks outside 
(inside) the magnetic Brillouin zone. Then, the narrow nearly localized band appears at the Hubbard gap 
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edge. This narrow band crosses the chemical potential for underdoped cases. Towards half filling, the velocity 
characterizing the narrow band appears to vanish. 

The Brinkman-Rice picture holds in that the quasiparticle mass diverges towards the MIT. As is well 
known, this is realized in the dynamical mean-field theoryES' which incorporates the Brinkman-Rice and 
Mott-Hubbard pictures. The theory gives nearly localized single-particle excitations at the chemical potential 
for doped cases. Since the momentum dependence of the self-energy is neglected, the dispersion of the 
Hubbard band takes basically simple one given by the Hubbard I approximation^ Then, it is natural that 
the resonant peak appears first near the (tt, tt) level for one-hole doping. However, this prediction does not 
reproduce the divergence of the charge compressibility. This is because in the dynamical mean-field theory, 
the peak of the local density of states shifts in the Hubbard band as the doping concentration changes. 
In contrast, as already mentioned, the sharp peak at the gap edge exists already in the insulating phase 
in our results and it pins the chemical potential. In addition, this one-hole state around (tt, tt) in the 
dynamical mean-field theory clearly violates the Luttinger sum rule. These problems are solved after AF 
spin correlations are properly included through the momentum dependence of the self-energy part. As we 
have already noted, AF spin correlations suppresses the dispersion outside (inside) the magnetic Brillouin 
zone in the lower (upper) Hubbard band and gives rise to formation of the SDW-like dispersions. 

The Slater picture holds only in the following two aspects; one, the ground state is antiferromagnetically 
ordered and the other, the SDW-like dispersions appear at low enough temperatures. This picture is improved 
when in addition to AF spin correlations, the growth of the local spin moment and thereby the splitting of 
the Hubbard bands are taken into considerations in the paramagnetic phase. This enables one to describe 
insulating state without the AF long-ranged order at finite temperatures as well as doped metallic states. 
Under strong AF spin correlations, the high-energy behaviour of the self-energy contributes to formation of 
the upper and lower Hubbard bands with the SDW-like dispersions. However, an extra narrow band appears 
at low energies as resonant single-particle excitations, in contrast with the Slater picture. 

§5. Comparison with high temperature superconductors 

Lastly, we discuss a strongly weak dispersion relation both around the (tt, 0) and (0, tt) momenta 
and along (tt, 0) Ux lar, tt) observed with the angle-resolved photoemission spectroscopy (ARPES) in 
Bi2Sr2CaCu208+5£j'E2P In the high-Tc compound, the dispersion around the (tt, 0) and (0,7r) momenta 
is significantly reduced compared with the band structure calculations. This flatness depends only weakly 
on both the temperature and the hole concentration; it persists below and above the superconducting critical 
temperature Tc- This suppression of the dispersion is attributed to strong correlation effects related to the 
MIT, as discussed previously based on QMC calculations .ES'LiP The present calculations of the 0PM gives 
an intuitive explanation of such enhancement of the electron mass: Local AF or singlet correlations increase 
|ia;,a;'| between the nearest-neighbor sets of x and x'. This affects the dispersion in the level of the two- 
pole approximation so that the (tt, tt) level in the lower Hubbard band shifts downwards and the dispersion 
around the (tt, 0) and (0, tt) momenta becomes flatter. Such strong correlation effects on the single-particle 
dispersions in the Hubbard bands can not be obtainciLif one introduces a decoupling approximation to the 
self-energy part as done in the FLEX approximation.cfP The similarity of the dispersions between the two- 
dimensional Hubbard model and high-Tc cuprates above the pseudogap temperature indicates that cuprates 
are actually located in the proximity of the MIT described by the Hubbard model. It also suggests that 
to capture essential features of the superconductivity, it is crucial to describe the correlation effects with 
strongly momentum-dependent renormalization obtained in this paper. We speculate that the low-energy 
narrow band obtained in this paper play an important role in realizing the superconductivity. Further studies 
are required to clarify its contribution to the superconductivity. 

§6. Summary 

We have developed the operator projection method (0PM) to study strong electron correlations in the Hub- 
bard model. The 0PM is a non-perturbative analytic-theoretical framework based on the continued-fraction 
or Dyson-equation expression of response functions.tiP It is remarkable that the method systematically im- 
proves the Hartree-Fock theory, approximations to the self-energy part based on the perturbation expansion 
or decoupling approximations, Hubbard approximations and two-pole approximations. In fact, the present 
theory captures three crucial ingredients of the MIT derived from (i) Mott-Hubbard, (ii) Brinkman-Rice and 
(iii) Slater pictures. Furthermore, our present theory incorporates the strong momentum dependence of the 
single-particle spectra, which emerges as a crucial aspect of the strong correlation and a source of electron 
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differentiation. 

This method reproduces the Mott-insulating phase as weU as the fiUing-control metal-insulator transition 
(MIT) in the Hubbard model in both one and two dimensions. A transition from metals to a Mott insulator 
occurs at half filling as a function of the filling. We have obtained the following characteristic features in the 
vicinity of the MIT in two dimensions: 

At half filling, the Mott gap separates two Hubbard bands. At low temperatures, the growth of AF 
spin correlations produces the AF shadows of the simple Hubbard bands. The formation of the Hubbard 
bands and their AF shadows result in a four-band structures at half filling. This four-band structure can be 
regarded as the superposition of two low-energy narrow bands and two SDW-like bands. The sharp peaks 
of the density of states at the gap edges are associated with an unusually flattened dispersion around (tt, 0) 
and (0, tt). At moderate temperatures, the AF shadows disappear and the single-particle excitations are well 
described by the Hubbard I approximation.B' At much higher temperatures T > C//2, incoherent low-energy 
spectral weights appear in the Mott gap, reflecting the loss of the quantum-mechanical features of electrons. 

In the hole-doped case, a narrow band appears at the top of the lower Hubbard band energetically separated 
from the upper Hubbard band by the Mott-Hubbard pseudogap. At low temperatures, a sharp peak of the 
local density of states develops near the Fermi level due to the narrow band. This low-energy narrow band 
has a particularly flat dispersion and dominant spectral weights around the (tt, 0) and (0, tt) momenta. This 
low-energy band is pinned at the edge of the incoherent lower Hubbard bands and generates low-energy 
hierarchy structure. Such extremely weak single-particle dispersion obtained in the present results agrees 
with the results observed in the high-Tc cuprates with the angle-resolved photoemission spectroscopy. The 
velocities characterizing the low-energy single-particle excitations vanish towards the MIT. This suggests 
that the quasiparticle mass diverges towards the MIT. The diverging compressibility is generated from the 
mass divergence of the pinned low-energy band at the Hubbard gap edge. This novel feature is certainly 
beyond the Brinkman-Rice picture. The damping rates for these low-energy excitations are also enhanced 
near half filling. Then, quasiparticles become ill-defined at least at T = 0.02t, since the damping rate remains 
comparable to wp- 

In the one-dimensional Hubbard model at half filling, the results for the single-particle spectra show similar 
Mott-insulating features for a Mott gap separating two Hubbard bands. As in two dimensions, strong AF spin 
correlations at low enough temperature produce AF shadows of the Hubbard bands in the plane of energy 
and momentum, in spite of the absence of the AF long-ranged order even at T = 0. Then, single-particle 
spectra have two branches in each Hubbard band. The bandwidth of the low-energy branch is narrower than 
the high-energy branch. This low-energy narrow branch is, however, not described as a simple spinon branch 
in the U = oo case. This is because the low-energy branch has finite weights for the momenta |A:| > tt/2 in 
the lower Hubbard band in the present results in contrast with the U — oo case. A mixing between spinous 
and holons is considered to yield such deviation from the U — oo results. 

The piffle main features obtained in the present study well reproduce the results of the QMC calcu- 
lations .HIIj^ For further improvements of the present 0PM, it is required to impose constraints on the 
consistency between the single-particle and two-particle properties. Furthermore, the resulting spin degrees 
of freedom at half filling in the strong-coupling regime must agree with predictions for the Heisenberg model. 
The higher-order projections are also important to extract low-energy properties more clearly. Studies on 
these improvements are currently in progress. 
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Fig. 2. The temperature dependence of the local density of states p{u!) for U/t = 1 (a), 4 (b) and 8 (c) at half filling (n) = 1. 
For (a), T/t = 0.02, 0.2, 0.5, 1 and 4. For (b), T/t = 0.02, 0.2, 0.5, 1, 4 and 8. For (c), T/t = 0.02, 0.5, 2, 5 and 8. At lower 
temperatures, less spectral weights remain around a; = 0. The Mott insulating gap grows at low temperatures T = 0.02t, 
while spectral weights gradually fill in the gap at high temperatures. It is remarkable that the Mott-Hubbard pseudogap 
remains even at T = J = 4t^ /U above which AFM spin correlations are not enhanced. On the other hand, at T > (7/2, a 
low-energy incoherent band develops. Detailed discussions are given in the text. 
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Fig. 3. Momentim| distribution»-^or C//t = 1 (a), 4 (b) and C//i = 8 (c) at half filling (n) = 1 compared with the QMC 
results from ref.E2P (b) and ref.H' (c). For comparison, momentum distributions for C/ = and (n) = 1 at the same set of 
temperatures are plotted with solid lines in the corresponding colors in each panel. The results support that the system is 
a non-Fermi liquid in all these cases. Although Fig. 2 shows that spectral weights develop around a; = at T > C//2, these 
low-energy excitations are too strongly damped to produce a jump in n^ . 
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Fig. 4. Single-particle spectral functions A(ijJ, k) = —Iin G{iLi, k)/n along (0, 0)— (tt, 0)— (tt, tt)— (0, 0) for U/t = 4 at half filling 
(n) = 1. The temperatures are T/t = 0.02, 0.5 and 4 in (a), (b) and (c), respectively. 
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Fig. 5. Single-particle dispersions along (0, 0)-(7r, 0)-(7r, tt)— (0, 0) for V jt = 4 at (n) = 1. The temperature is T/i = ftT02, 0.5 
and 4 in (a), (b) and (c), respectively. The squares denote the present results and the red circles denote QMC results.|y The 
open and the filled symbols represent the momenta where the peak intensities divided by the largest intensity are less mA 
larger than 0.1, respectively. Main features of our results marked by the filled symbols are similar to those of QMC resultsjf^ 
except that the Mott-insulating gap is overestimated around the (tt, 0) and (0, tt) momenta in the present results. The solid 
lines denote the free-electron dispersions. The dotted and the dashed lines represent the dispersion of the upper Hubbard 
band within the Hubbard I approximation and that obtained by the two-pole approximations with t^ = t^.Q. 
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Fig. 6. The temperature dependence of the local density of states p{lo) for U/t = 1 (a), 4 (b) and 8 (c) at n = 0.87. For (b), 
T/t = 0.02, 0.2, 0.5, 1 and 4. For (c), T/t = 0.02, 0.5, 2, 5 and 8. A sharp peak develops near o) = at low T. 
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Fig. 7. Momentum distributions along (0, 0)— (tt, 0)-(7r, 7r)-(0, 0) for C//t = 1 (a), 4 (b) and 8 (c) at (n) = 0.87. For comparison, 
momentum distributions for (7 = and (n) = 0.87 at the same set of temperatures are plotted with solid lines in the 
corresponding colors in each panel. The Fermi momenta fcp where n^ most strongly changes move towards (tt, tt) with 
increasing temperature. Thus, the "Fermi volume" increases with increasing temperature. 
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Fig. 8. Single-particle spectral functions A{uj, k) = -Im G{uj, k)/iT along (0, 0)-(7r, 0)-(7r, 7r)-(0, 0) for U/t = 8 at {n> = 0.87. 
The temperature is T/t = 0.02, 0.5 and 5 in (a), (b) and (c), respectively. They show that around the (7r,7r) momentum, the 
dispersion of the lower Hubbard band becomes weak and the weights also become small. The dispersion around the (tt, 0) 
and (0, vr) momenta are also found to be suppressed. 
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Fig. 9. Single-particle dispersions along (0, 0)-{7r, 0)-{7r, 7r)-{0, 0) for U/t = 8 at (ra) = 0.87. The temperature is T/t = 0.02, 
0.5 and 5 in (a), (b) and (c), respectively. The open and the filled symbols represent the momenta where the peak intensities 
divided by the largest value of A(ui,k) arc less and larger than 0.1, respectively. The solid lines denote the free-electron 
dispersions. 
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Fig. 10. Single-particle spectral functions A{uj, k) = -Im G(w, fc)/7r along (0, 0)-(7r, 0)-(7r, 7r)-(0, 0) for U/t = 4 at (n) = 0.87. 
The temperature is T/t = 0.02, 0.5 and 4 in (a), (b) and (c), respectively. 
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Fig. 11. The chemical potential measured from the particle-hole symmetric energy /i — U/2 versus the hole concentration 
1 — (n) for the two-dimensional Hubbard model with U/t = 1,4 and 8 at T/t = 0.02. 
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Fig. 12. Fermi velocities vp and damping rates 7^ are given as a function of the hole concentration 1 — (n). They are given 
for the Fermi momenta in directions along both (0, 0)-{7r, 0)-{7r, vr) and (0, 0)-(7r, tt). 7j^ is determined as the width of the 
lorentzian fit to A(lj, k) at T/t = 0.02. The figure shows that in both momentum directions, the quasiparticle mass diverges 
towards the MIT, as in the Brinkman-Rice picture. However, v/^ exhibits a momentum dependence so that that around 
(vr, 0) and (0, tt) is less than that along (0, 0)-(7r, tt). It also shows the enhancement of the quasiparticle damping rates as 
(n) — » 1. The lines in this figure are guide to the eye. In the case oi U = 0, vp/t = 2 and 1.998 along the momentum 
direction (0, 0)-(7r, tt) at (n) = 1 and 0.87, respectively. Along (0, 0)-(7r, 0)-(7r, tt), vp/t = and 0.942 at (n) = 1 and 0.87. 
For 0.94 < (n) < 1, low-energy excitations near the Fermi level and momenta in A{ijj,k) exhibit a too broad structure to 
recognize quasiparticles with well-defined velocities and damping rates. Therefore, the Fermi surface is also ill-defined in this 
filling region at least at this temperature T/t = 0.02. 
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Fig. 13. Spectral function A{ll), k) as a function of the momentum in the one-dimensional Hubbard model for U /t ■ 
filling. The temperature is T/t = 0.01 (a) and 0.2 (b). 
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